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Abstract. We solve the unstructured search problem in constant time by computing 
with a physically motivated nonlinearity of the Gross-Pitaevskii type. This speedup 
comes, however, at the novel expense of increasing the time-measurement precision. 
Jointly optimizing these resource requirements results in an overall scaling of TV 1 / 4 . 
This is a significant, but not unreasonable, improvement over the N 1 / 2 scaling of 
Grover's algorithm. Since the Gross-Pitaevskii equation approximates the multi- 
particle (linear) Schrodinger equation, for which Grover's algorithm is optimal, our 
result leads to a quantum information-theoretic lower bound on the number of particles 
needed for this approximation to hold, asymptotically 

PACS numbers: 03.67.Ac, 05.45.-a, 67.85.Hj, 67.85. Jk 
1. Introduction 

Abrams and Lloyd pQ argued that a nonlinear quantum theory could result in 
unreasonable computational advantages by giving two examples of nonlinear algorithms 
that solve NP-complete and #P problems in polynomial time. Both of their algorithms 
can be implemented by a nonlinear Schrodinger-type evolution in which the time 
derivatives of the state components depend upon their hyperbolic tangents [21 Ej- The 
derivative of tanhx at x — is 1, so this is a strongly nonlinear system in which is an 
unstable fixed point. The strength of the nonlinearity provides a large computational 
advantage, but it also makes the system highly susceptible to noise [U 121 E]- 

An obvious question is whether a modest, physically motivated nonlinearity can still 
produce a computational advantage. In particular, consider Bose-Einstein condensates. 
In 1924-25, Bose and Einstein predicted that cooling a dilute gas of bosons near absolute 
zero would cause the atoms to occupy their lowest quantum state, forming a new 
state of matter where quantum effects would be macroscopically apparent (H EJ Ej- It 
took seventy years, however, for these Bose-Einstein condensates to be experimentally 
produced [3 El E] • In general, describing such many-body systems is difficult because 
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of the many interaction terms. But under certain conditions, one can assume that only 
two-body contact interactions contribute and the s-wave scattering length a is much 
less than the interparticle spacing. Then using mean field theory, one finds that the 
system is approximately described by a nonlinear Schrodinger equation with a cubic 
nonlinearity: 

^M), (i) 

where H includes the kinetic energy and trapping potential, m is the mass of the 
condensate atom, and N is the number of condensate atoms||] The validity of this 
celebrated Gross-Pitaevskii equation [TOf [TT] requires that N be much greater than 
1 — but how much greater? 

We provide a quantum information-theoretic solution to this question in the context 
of solving the unstructured quantum search problem [T2] in continuous-time [32] using 
the Gross-Pitaevskii equation as the governing equation. The cubic nonlinearity in ((!]) 
has zero derivative at zero, making it softer than those considered by Abrams and Lloyd 
PQ . We quantify the computational advantage that such a nonlinearity provides for the 
unstructured search problem compared to standard quantum computation; this requires 
considering time-measurement precision as a physical resource. Since this advantage 
cannot persist when the Gross-Pitaevskii equation is recognized as an approximation 
to an underlying multi-particle Schrodinger equation, for which Grover's algorithm is 
optimal, we arrive at a quantum information-theoretic lower bound on the number of 
particles, N , needed for this approximation to hold, asymptotically. 

2. Setup 

The system evolves in a ^-dimensional Hilbert space with computational basis 
{|0),...,|iV — 1)}. The initial state \ip(0)) is an equal superposition \s) of all these 
basis states: 

N-l 

l*<o» = W = ^X> 

The goal is to find a particular "marked" basis state, which we label \w). 

Let's first review the linear solution. We use Childs and Goldstone's [2] notation 
for Farhi and Gutmann's [13] Hamiltonian: 

H = — jN\s)(s\ — 

where 7 is a parameter, inversely proportional to mass. The system evolves in the two- 
dimensional subspace spanned by \w) and \s). One might (correctly) reason that the 
success of the algorithm depends on the value of 7. This can be seen in figure [TJ which 
shows the difference in eigenvalues of — Ho and the overlaps of its nontrivial eigenvectors 
with \s) and \w). When 7 takes a critical value of 7 C = 1/N, the eigenstates of — H 

I More generally, the cubic nonlinear Schrodinger equation is the equation of motion for 4 field theory. 
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Figure 1. Eigenvalue gap and eigenstate overlaps of — Hq with TV = 1024. 



are proportional to ±|w) + \s), and the corresponding energy gap is 2/yN. So the 
Schrodinger evolution rotates the state from \s) to \w) in time tt\/N /2. 

In the nonlinear regime, we include an additional nonlinear "self-potential" V(t) so 
that the system evolves according to the Gross-Pitaevskii equation Q: 

i^(r,t) = [H - g\i>(r,t)\ 2 W(r,t), 
at v ' 

V(t) 

where g > 0. This corresponds to a BEC with attractive interactions, and thus a 
negative scattering length [Tj)J US]- Heuristically, as probability accumulates at the 
marked state due to the term in H , the self-potential attracts more probability, 

speeding up the search. Thus we expect larger g to result in a faster algorithm. 
In the computational basis, the self-potential is 

N-l 

v{t)= g 52\m\ 2 \m 

i=0 

Even with this nonlinearity, the system remains in the subspace spanned by {\w), \s)} 
throughout its evolution. We define a vector 



|r) 



which is orthonormal to \w). Then the state of the system \ip(t)) can be written as 

m)) = a(t)\w)+/3(t)\r). 
Writing the Gross-Pitaevskii equation in this {\w), |r)} basis, we get 

d 



a 

dt \ (3 



-i Hn 



V) 



a 





7 + 1 + g\a\ 

^y/N - 1 



7(iV-l) + ^|/3| 2 ) { (3 



(2) 



.4 
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Figure 2. Success probability as a function of time for N = 1024 and 7 = 1/N 
constant. The solid curve is the linear (g = 0) case, and the dashed curve is the 
nonlinear ,g = l case. 




Figure 3. Eigenvalue gap and eigenstate overlaps of A at t = 20 for nonlinear search 
with N = 1024, g = 1, and 7 = 1/N constant. 



where we've defined A = —(H — V). 
3. Critical Gamma 

Before proceeding with further analytical calculations, we build some intuition by 
examining two plots. For constant 7 and g, the success probability as a function of 
time, |a(t)| 2 , is plotted in figure [2] along with the linear result. The nonlinear algorithm 
underperforms the linear one in this case. This is true in general for constant 7 and g, 
and it can be understood by examining the time-dependence of the critical value of 7, 
which is the value of 7 that ensures that the eigenstates of A are in the form ±|u>) + \ s). 
Initially, 7 C = 1/N. Then, as shown in figure [3j it shifts to a larger value. If 7 is 
constant, it will not follow this shift, we will no longer have the desired eigenstates, and 
the algorithm will perform poorly. 
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To determine how 7 C varies with time, we find the eigenvectors of A and choose 7 
so that they have the desired form ±\w) + \ s) . To eliminate fractions in the subsequent 
algebra, we rescale the nonlinearity coefficient g by defining 

Solving the characteristic equation gives the eigenvalues of A: 



A ± = -( 7 iV 



l + Ga)±^AA, 



where the gap between them is 

AA = y/frN - l) 2 + 4 7 + G 2 5 2 + 2G5 [1 - 7 (iV 
and we've defined 



2)], 



a = (N- l)\a\ 2 + \(3\ 2 and 5 = (N — l)\a\ 2 — \f3\ 2 . 



The corresponding eigenvectors of A are 

N f-7iV + 1 + 5G± AA 



N 



w) + \s) 



2^VN 

The critical value of 7 ensures that these eigenvectors have the form ±|if) + \s). That 
is, 

-7JV + 1 + 5G ± AA 



Solving this yields: 



7c 



G5 



1. 



7c 



N 



(3) 

Note that in the linear limit (G = 0), this reduces to 7 C = 1/N, as expected. 
Importantly, since S varies with time, ^ implies 7 C also varies with time, in agreement 
with our previous discussion about figures [2] and [3j 



4. Runtime 

For the remainder of the paper, we choose time- varying 7 = j c according to (|3]). Before 
analytically determining the consequences of this, let's again consider a plot. Figure [4] 
shows the success probability as a function of time. There are several observations. First, 
the success probability reaches 1, which occurs because we constructed the eigenstates 
to make this happen. Second, as N increases, the runtime remains constant. Third, 
the success probability is periodic. Finally, the peak in success probability becomes 
increasingly narrow for large N. Let's now analytically prove the second, third, and 
fourth observations. 

To begin, we explicitly write out ^ to get two coupled, first-order ODEs for a(t) 
and /3 (t): 

^ = i { [7c + 1 + G(N - l)\a\ 2 } a + lc VN^l^ (4) 
^ = i { lc VN^la + [ lc (N - 1) + G|/3| 2 ] /?} . (5) 
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Figure 4. Success probability as a function of time for nonlinear search with G = 1 
and 7 = 7 C as defined in The solid line is N = 100 and the dashed line is N — 1000. 



We can decouple these equations for the success probability x(t) = |a(t)| 2 , yielding 



dx 
~dt 



U{Nx - 1)(1 - x) [1 + G(iVx - l)f 



N 2 



(6) 



The details of this decoupling procedure are given in |Appendix A[ To solve this 
uncoupled equation, we use separation of variables and integrate from t = to t and 
x = 1/N to x, which yields 



t 




N 

G(N-V 



tan 



-i 



y/l + G(N-l)y/N. 



X 



1 



7T 
2 



(7) 



Solving for x, the success probability as a function of time is 
x{t) 



iV+[l + G'(iV-l)]tan 2 


/l +G(A r_i) 
2 V N L 


iV + N [1 + G(N - 1)] tan 2 


/l+G(JV-l) . 
2 V 



(8) 



From this, the success probability reaches 1 when the tangent term is zero, which first 
occurs at time 

1 tiVN 



v/l + G(iV-l) 2 

This runtime is exactly constant for G = 1. Also, when G = 0(1), the runtime for large 



is 7r/2vG, and thus asymptotically constant (and arbitrarily small!). From (|8j), we 
also see that the success probability is periodic with a period of 2t*. 

Now let's prove that the peak in success probability is narrow by finding its width, 
thus proving all our observations about figure |1J Using ([7]), the difference in time at 
which the success probability reaches a height of 1 — e is 



At 



N 



1 + G(N-1] 



tan 



^l + G(N-l)y/N{l-e)-l 
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The tan 1 makes it difficult to determine the scaling with N, so we Taylor expand it: 



^ITcfc)V& 0(f3/2) ' 

When G = N K , the first term scales as 9(A^ 1 / 2 ) when k < — 1 and 9(A r ~ 1 / 2_K ) when 
K > — 1, for large N. To determine whether keeping this first term alone is sufficient, 
we use Taylor's remainder theorem to bound the error 



N*(l + 3G(N(l-e)-l)) ^ 
(N(l - e) - if 2 (1 + G (N (1 - e) - l)) 2 



which has the same scaling for large N as the first term in the Taylor series for At. 
Thus it suffices to keep only the first term. 

For constant G, the width in success probability is 0(l/vW), which agrees with 
our observation from figure [4] that the peak in success probability is increasingly narrow 
as N increases. Thus we must measure the system with increasing time precision. This 
behavior is opposite the linear case. That is, when G = the width is O(VN), so the 
time at which we measure the result can be increasingly imprecise as N increases. 



5. Time-Measurement Precision 



This time-measurement precision requirement of the nonlinear algorithm requires 
additional resources. In particular, time and frequency standards are currently defined 
by atomic clocks, such as NIST-F1 in the United States [17]. An atomic clock with 
^ciock ions used as atomic oscillators has a time-measurement precision of 1 / ^/^ciock 
when the ions are acted upon independently. This can be improved using quantum 
entanglement, reducing the time-measurement precision to 1/ridock [18, 19J. Even with 
this improvement, our constant-time nonlinear search algorithm would require 0(\/N~) 
ions in an atomic clock to have sufficiently high time-measurement precision to measure 
the peak in success probability. So, although our nonlinear algorithm runs in constant 
time, the total resource requirement is still 0(\/iV), the same as the linear algorithm. 
This raises the possibility that nonlinear quantum mechanics may not provide efficient 
solutions to NP-complete and #P-complete problems when all the resource requirements 
are taken into consideration pQ. 

In our case, however, we can settle for a smaller improvement in runtime and reduce 
the time-measurement precision and total resource requirement. If we let G decrease as 
N K for k < 0, then the runtime is t* = 6(iV~ K / 2 ), and the time-measurement precision is 
At = 0(iV -1 / 2-K ), where we've assumed for both that k > —1, since for k < —1, At = 
0(iW 2 ), independently of G. This time-measurement precision requires 0(A^ 1 / 2+K ) ions 
in an atomic clock that utilizes entanglement. We assume, as in the setup for Grover's 
algorithm, that logiV qubits can be used to encode the N- dimensional Hilbert space; 
these should also be included in the required "space" resources. Multiplying the time 
and "space" requirements together, which preserves the time-space tradeoff inherent in 
naive parallelization, the resulting total resource requirement takes a minimum value of 
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Figure 5. Success probability as a function of time for nonlinear search with 
G = N~ 1/2 and 7 = 7c as defined in The solid line is N = 100 and the dashed 
line is N = 1000. The peaks have same width, independent of N. 



0(N 1 ^ log N) when n = —1/2 (so that the runtime is iV 1 / 4 and the time-measurement 
precision is constant). The success probability as a function of time at this jointly 
optimized value of G is plotted in figure [5j note that the peak width is independent of 
N. 

This significant — but not unreasonable — improvement over the Q(^/N log N) time- 
space resource requirements of the linear quantum search algorithm is consistent with 
our expectation that a modest nonlinearity should result in a modest speedup. 

6. Repulsive Interactions 

To prove the periodicity of the success probability x(t), we set ^ equal to zero and 
find critical points x* = 1/N, 1, and (G — 1)/NG. Then by examining the second 
derivative of x(t) at these critical points (either by differentiating ^ or using (A.9) 



that naturally arises in the decoupling procedure), we find that the critical points are 
minima, maxima, and stationary points, respectively. The stationary point is between 
1/N and 1 when G < —1/(N — 1), so it is irrelevant for attractive interactions since 
G > 0. Thus x(t) begins at the minimum 1/N and takes the positive square root in ^ 
up to the maximum value of 1. Since ^ is autonomous, meaning it only depends on 
x and not t, the success probability then takes the negative square root down until it 
again reaches the minimum 1/N. This cycle repeats, so x(t) is periodic. 

Our nonlinear search algorithm was based on the intuition that attractive 
interactions speed up the accumulation of success probability. By the same intuition, 
repulsive interactions, where G < 0, should yield a worse runtime. Our derivation of 
(|3j) for the critical value of 7 is unchanged if we flip the sign of G, so JO) and @ are 



still valid for repulsive interactions. As previously stated, these equations yield critical 
points x* = 1/N, 1, and (G — 1)/NG, corresponding to minima, maxima, and stationary 
points, respectively. 
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Figure 6. Success probability as a function of time for TV = 1024 and 7 = 7 C as 
denned in The solid line is the linear (g — 0) case, the dashed line is the nonlinear 
g = —0.5 case, the dotted line is the nonlinear g = — 1 case, and the dot-dashed line is 
the nonlinear g = —1.5 case. 

When G > — l/(N — 1), the success probability is unhindered by the stationary 
point and reaches a maximum value of 1, as shown in the dashed curve of figure [6j 
When G < — 1/(N — 1), however, reaching this maximum is precluded by the presence 
of a stationary point, as shown in the dashed and dot-dashed curves of figure [6] 

We can explicitly prove that repulsive interations (G < 0) will underperform the 
linear (G = 0) algorithm. From ([6]), 

^ = ±^V(Nx-l)(l-x)[l + G(Nx - 1)] . 

So when G < 0, the magnitude of dx/dt at a particular value of x is less than when 
G = 0. Then success probability will increase more slowly for repulsive interactions 
(except initially, where they increase at the same rate). Thus it will underperform the 
linear algorithm. 

7. Validity of the Gross-Pitaevskii Equation 

Of course, the cubic nonlinearity we've exploited is not fundamental, but rather occurs 
in an effective description of an interacting multi-particle quantum system (e.g., a 
BEC). So we must include the number of particles N Q in our resource accounting. 
Each particle interacts with the potential at the marked site, so in the framework of 
Zalka's optimality proof for Grover's algorithm [20] (generalized to continuous time 
[21]), there are N oracles, each responding to a logiV bit query. Zalka showed that 
the product of the space requirements and the square of the time requirements is 
lower bounded by N, i.e., (N logiV)(A^ 1 / 4 ) 2 = Cl(N). Solving for the number of 
particles, iVo = ^A^/logiV). This 

is a quantum information-theoretic lower bound 
on the number of particles necessary for the Gross-Pitaevskii equation to be the correct 
asymptotic description of the multi-particle (linear) quantum dynamics. 
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Notice that once we account for the scaling of Nq in the space requirements, the 
product of the time and space requirements is 0(iV 3 / 4 ), worse than the 0(N^ 2 log N) 
of Grover's algorithm. In fact, if we calculate for the general case G = N K , where 
k need not be chosen to optimize the product of the time and space (ignoring N ) 
resources, Zalka's bound implies N = fi(max{l, N 1+K / log N}), so the total time-space 
requirements are 0(N 1+K ^ 2 ) for k > — 1, and 0(N 1 ^ 2 log N) when k = —1. This is 
optimized for k = —1, i.e., by Grover's algorithm. On the other hand, Zalka's bound 
is strongest when k = 0, in which case it implies that Nq = Q(N/\ogN). That is, the 
existence of the constant time nonlinear algorithm we found in section 4 implies this 
stronger lower bound on Nq, despite the 0(N 1 ^ 2 ) number of clock ions required. To our 
knowledge, this is the first lower bound derived on the scaling of Nq required for the 
Gross-Pitaevskii equation be a good asymptotic approximation. 

This bound also is significantly stronger than the bound implied by the physically 
plausible requirement that the volume of the multi-particle condensate, and thus N , be 
of at least the order of the volume of space in which the N possible discrete locations 
are defined. Were we working in any fixed, finite dimension, e.g., on a cubic lattice, 
the volume would be proportional to N, implying N = Q(N). But we are not; the 
complete graph with equal pairwise transition rates is realized by the vertices and edges 
of an equilateral (N — l)-dimensional simplex. With edges of length 1, this has volume 
y / N/2 N ~ 1 /(N — 1)!, which is much smaller than N, and also much smaller than our 
bound of N / log N. 



8. Critical Gamma is a Continuous Rescaling of Time 

We previously derived the critical value of 7 so that the eigenstates of the Hamiltonian 
are proportional to ±|iw) + \s). Now we examine what the critical value of 7 does from 
another perspective. Recall the "Hamiltonian" we've been using is 

H = - 1 N\s)(s\ - \w)(w\- gJ2\H 2 \i)(*l 

i 

where ipi = Explicitly writing the nonlinear term as marked and unmarked terms, 

we get 

H = -jN\s)(s\ - \w)(w\ - g\a\ 2 \w)(w\ - g — - ^ \x)(x\ 

= —rN\s)(s\ - \w)(w\ - G(N - l)\a\ 2 \w)(w\ - G\f3\ 2 ^ \x}(x\ 
= -yN\s)(s\ - [1 + G(N - l)\a\ 2 } \w)(w\ - G\f3\ 2 ^ \x)(x\. 

x^w 

Recall 7 = 7 C is chosen according to (j3|: 

7c iV = l + G(iV-l)|a| 2 -G|/3| 2 , 
which we arrange to get 

1 + G(iV - 1)H 2 = 7c iV + G|/3| 2 . 
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Then the Hamiltonian becomes 

H = -j c N\s)(s\ - [ lc N + G\f3\ 2 ] \w)(w\ - G\f3\ 2 ^ \x)(x\ 

= - lc N(\s)(s\ + \w)(w\)-G\P\ 2 I. 

The last term continuously redefines the "zero" of energy, so we can drop it. That is, 
it only changes the overall phase of the system, which has no measurable effect. Then 
the Hamiltonian is 

H =-^N(\s)(s\ + \w)(w\). 

Importantly, Hfq = —\s)(s\ — is the Hamiltonian from Farhi and Gutmann's 

"analog analogue" of Grover's algorithm [13], and it is optimal. Our nonlinear algorithm 
has a factor of jN, so it effectively follows their optimal algorithm, but with a 
continuously rescaled time. That is, the system evolves according to 

Let's call the rescaled time r(£) so that dr = ^Ndt. Then 
7iVdt, 



and the equation of motion becomes 

1— = H FG ip. 
dr 

This has success probability given by (11) of [13J: 

. ) / / \ 1 o ( t~ 

x[t) = sin — — + — cos 



N N \Jn 



Plugging in for r, 



Since 7 C A^ = 1 + GS=1-G + GNx, we get 



x{t) = sin 



1 - G)t + GN Jx(t)dt 



N 



1 „ /(1-G)t + GN [x(t)dt 
+ — cos 1 



This integral transcendental equation gives x(t). While the difficulty of solving this 
equation makes it less useful in practice, it does reveal our nonlinear algorithm's 
relationship with the linear, optimal algorithm. In particular, a different control policy 
for 7 will cause the system to evolve along a different, slower path. While not a proof, 
this is an argument for the optimality of our algorithm. 
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9. Multiple Marked States 



Our analysis naturally extends to the case of k marked states. Let M be the set of 
marked basis states. As before, the system evolves in a two-dimensional subspace: 



* x<£M 



The system evolves according to 
d 



dt 



\m) = iA\ 



where 



A = jN\s)(s\ + 



xeM " x£M 

includes both the linear Hamiltonian and the nonlinear "self-potential" . The eigenstates 
of A have the form ±|u>) + \s) when 7 is 



7c 



1 + G6 



where G = g/(k(N — k)) and 5 = (N — A;)|a| 2 — k\(3\ 2 . At 7 = j c , we can decouple these 
equations in the same manner as the k — 1 case (see Appendix A) and integrate from 
t = to t and x = k/N to x to get 



t 



N 



tan 



-1 



k(l + G(N — k)) 
which can be solved for a success probability of 

x(t) 



a/1 + G(N - Fjy/Nx - k 



TV 

2 



N + k[l + G(N - 


k)] tan 2 


7T 

2 


^Jk(l+G(N-k)) t 


N + N [1 + G(N - 


- k)] tan 2 


7T 

2 


/k(l+G(N-k)) t 



Then the runtime is 



y/k{l + G{N-k)) 2 ' 
and the success probability is still periodic with period It*. At this runtime, the peak 
in success probability has a width of 



At 



N 



tan 



'l + G(N-k)y/N(l-e)-k 



k(l + G(N - k)) 

but Taylor's theorem can be used to show that it suffices to keep the first term in the 
Taylor series: 

2N 



At 



+ 0(e 



3/2\ 



l + G(N-k) V k(N-k) 

As in the case of a single marked state, we can find the scaling of G = N K that 
optimizes the product of "space" and time, where "space" includes both the number 
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of ions needed in an atomic clock that utilizes entanglement to achive sufficiently high 
time-measurement precision, and the log N qubits needed to encode the iV-dimensional 
Hilbert space. Say the number of marked sites scales as k = N x , with < A < 1. 
When k = —A/2 — 1/2, the product of "space" and time takes a minimum value of 
ST = jV~ A//4+1//4 log iV (so that the runtime is jV~ A / 4+1 / 4 and the time-measurement 
precision is constant). Note this is a square root speedup over the linear (G = 0) 
algorithm, whose product of "space" and time is jV~ A / 2+1 / 2 log N. Thus our nonlinear 
method, by varying 7 and choosing an optimal nonlinear coefficient G, provides a 
significant, but not unreasonable, improvement over the continuous-time analogue of 
Grover's algorithm, even with multiple marked items. 
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Appendix A. Decoupling the Equations of Motion 



We begin with two coupled, first-order ODEs 
da 



dt 

d/3 

dt 



i { [ 7c + 1 + G(N - l)\a\ 2 ] a + 1c VN~^1/3^ 
i { 7c v / iV^Ta + [ic(N - 1) + G\P\ 2 ] f3 



(A.l) 
(A.2) 



We decouple these equations by defining three real variables x(t), y(t), and z{t) such 
that 



x = \a\ 

y + iz = a/3*. 



(A.3) 
(A.4) 



Note that x{i) defined by (A.3) is the success probability. Differentiating it and utilizing 

da* 



(A.l ), we find that 

dx d\a\ 2 



dt dt 
Solving this for z, we get 

1 



a 



— h —a 

dt dt 



2 lc VN - lz. 



dx 



2 7cV / A r - 1 dt 



Noting that d7 c /dt = Gdx/dt, we differentiate (A. 5) to get 



dz _ 1 
dt ~ 2 v / iV^T 



1 „ { dx \ 2 1 d 2 x 



-G 



+ 



dt J 7 C dt 2 



(A.5) 



(A.6) 
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Now we want to find another expression for dz/dt, which we can then set equal to (A. 6). 
We do this by differentiating (A.4), utilizing (A.l) and (A.2), and equating the real and 
imaginary parts, which yields 
dy 

dt 

dz 

dt 



-2-7 c z 

2 7c y + 7c v / A- 1(1 -2x). 



(A.7) 
(A.8) 



Substituting (A. 5) for z into (A.7), we get 

dy — 1 dx 
dt ~ v^V^Idt' 

which integrates to 



1 



y 



where the constant of integration was found using x(0) = l/N and y(0) = y/N - l/N. 
Now we can plug this into (A.8) to get 
dz _ 1 
dt 



27c 



7c 



(1 - x) + 7 C VN- 1(1 - 2x) 



[1 + N- 2Nx) . 



Equating this to (A. 6) and simplifying yields 

d 2 x G ( dx\ 2 „ , , T „ r . 
dF = ^U) + ^(1 + ^-2^). 

Plugging in for 7 C as defined in (J3l), this becomes 



d 2 x 
dt 2 " 



NG 



dx\ 



1-G + NGx V dt J 
o the 
NG 



+ ^-(l-G + NGx) 2 (1 + N- 2Nx) (A.9) 
N 2 



Now let /(x) = (dx/dt) 2 so that d//dx = 2d 2 x/dt 2 . Then ( |AT9[ ) becomes 

Id/ 

2dx 



Solving this first-order ODE and using the initial condition f(x = l/N) = 0, we get 

4(Ax - 1)(1 - x) [1 + G(Ax - l)] 2 



/(*) 



A 2 



Taking the square root and noting that dx/dt = ±\J~fi 
± 



dx 
dt 



'4(Ax - 1)(1 - x) [1 + G(Ax - l)f 



A 2 



(A.10) 
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